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Disclaimer 


- These sheets contain mathematics to illustrate properties of Markov 
Random Fields (undirected networks) 


- As well as for reference 


- Without mathematical statistics some concepts might be hard to 
follow 


- And that is ok! 


- | will not ask you to repeat the mathematics The main thing you 
should take from this lecture is: 


- Understand what Markov Random Fields are 
- Understand what independence relations they imply 
- Understand how the edge weights should be interpreted 
- Understand how they are estimated 
- As we will see on Thursday, the R codes are simple!!! 
Be there on Thursday! 
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Recap 


Causality 


Easiness of Class Intelligence 
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Causal structures 


Chain Common Cause Collider 


Intelligence Easiness of Class 


Intelligence 
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Easiness of Class Intelligence 
e 


What if we don't know the structure? 
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Easiness of Class |-- Intelligence 
e 


Are the two nodes independent given any set of other nodes (including 
the empty set)? 


- Yes! They are independent to begin with! 
Draw no edge between Easiness of Class and Intelligence 
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Easiness of Class Intelligence 
‘“ 
‘“ 
‘“ 


Are the two nodes independent given any set of other nodes (including 
the empty set)? 


No! 
Draw an edge between Easiness of Class and Grade 
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Easiness of Class Intelligence 


Are the two nodes independent given any set of other nodes (including 
the empty set)? 


No! 
Draw an edge between Grade and Intelligence 
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Easiness of Class Intelligence 
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Easiness of Class Intelligence 






Is the middle node in the set that separated the other two nodes? 


- Yes! 
No nothing 
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Easiness of Class Intelligence 







Is the middle node in the set that separated the other two nodes? 


No! 
Grade is a collider between Easiness of Class and Intelligence 
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Easiness of Class Intelligence 







Do we now know the direction of the edge between Intelligence and 
IQ? 
- No! 
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Easiness of Class Intelligence 







Do we now know the direction of the edge between Grade and 
Diploma? 


- Yes! Grade was not a collider! 
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Easiness of Class Intelligence 
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Equivalent models 


Easiness of Class Intelligence Easiness of Class Intelligence 
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Outline 


- Testing Directed Networks 

- Structural Equation Modeling 

- Problems with directed networks 
Undirected networks / Markov Random Fields 
Undirected Networks as generating structures 
Important Markov Random Fields 

- Ising Model 

- Gaussian Random Field 
Examples 

- Radicalization 

- Quality of Life 
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Testing Directed Networks 


Determinants of Radicalization of Islamic Youth in the 
Netherlands: Personal Uncertainty, Perceived 


Injustice, and Perceived Group Threat 


Bertjan Doosje* 
University of Amsterdam 


Annemarie Loseman and Kees van den Bos 
Utrecht University 


In this study among Dutch Muslim youth (N = 131), we focus on the process of 

radicalization. We hypothesize that this process is driven by three main factors: (a) 

personal uncertainty, (b) perceived injustice, and (c) perceived group threat. Using 

structural equation modeling, we demonstrate that personal uncertainty, perceived 

injustice, and group-threat factors are important determinants of a radical belief 

system (e.g., perceived superiority of Muslims, perceived illegitimacy of Dutch 

authorities, perceived distance to others, and a feeling of being disconnected 

from society). This radical belief system in turn predicts attitudes toward violence 

by other Muslims, which is a determinant of own violent intentions. Results are 

discussed in terms of the role of individual and group-based determinants of 

radicalization. 
Doosje, B., Loseman, A., & Bos, K. (2013). Determinants of radicalization of Islamic youth in the 
Netherlands: Personal uncertainty, perceived injustice, and perceived group threat. Journal of 


Social Issues, 69(3), 586-604. 
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Table 1. The Means, Standard Deviations, and Inter-Correlations of All the Constructs 





M 
1. Identification 4.56 
2. Ind. Rel. Depri. 2.39 
3. Col. Rel. Depri. 3.31 
4. Int. Anxiety — 0.20 
5. Symbolic Threat 3.46 
6. Realistic Threat 3.10 
7. Per. Em. Uncertain. 2.84 
8. Perc. Proc. Injustice 2.38 
9. Perc. Illegitimacy 2.37 
10. Perc. Ingr. Super. 3.26 
11. Distance 2.32 
12. Disconnected 2.79 
13. Moslim Violence 2.89 
14. Violent Intentions 2.08 


SD 


0.85 
0.81 
0.92 
0.17 
0.76 
0.88 
0.67 
0.68 
0.02 
0.93 
0.66 
0.96 
1.06 
0.91 
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Note. 2 = Individual Relative Deprivation, 3 = Collective Relative Deprivation, 4 = Intergroup Anxiety, 5 = Symbolic Threat, 6 = Realistic Threat, 7 = 
Personal Emotional Uncertainty, 8 = Perceived Procedural Injustice, 9 = Perceived Illegitimacy, 10 = Perceived Ingroup Superiority. *p < .05. 
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Background Determinants 
Variables 





R?=17 
Intergroup 
Anxiety 
In-group 
Identification 


Individual 
Deprivation 


Collective Realistic 
Deprivation Threat 


R*=.07 
Personal 
Emotional 
Uncertainty 


R?=.27 
Perceived 
Injustice 


Radical Belief 
System 


_R=18 
Perceived 


In-group 
Superiority 


R°=.05 
Distance to 
Other People 


R?=.05 
Societal Dis- 
Connected 





Violent 
Intentions 


_ R2=39 
Attitude 
towards Muslim 
Violence 


R?=.26 
Own Violent 
Intentions 


Fig. 1. Final structural equation model. All paths are significant. R? = % variance explained. 
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The hypothesized model had a good fit: Chi-square (65) = 76.58, p = .154, 
CFI = .98, NFI = .87, GFI = .93, SRMR = .082, and RMSEA = .037. La Grange 
Multiplier Test suggested including two direct paths: from collective deprivation 
to perceived illegitimacy of Dutch authorities, and from perceived distance to own 
violent intentions. When we included these paths, the fit became better. Our final 
model is presented in Figure 1. It has a very good fit: Chi-square (62) = 58.13, 
p = .650, CFI = 1.00, NFI = .90, GFI = .94, SRMR = .070, and RMSEA = .000. 
All paths included in the model are significant. We discuss this model in steps 
from left to right. 


25/110 


Equivalent model 


vs 
———— — 
col_dep|— 





| | \ i | 

\ injustice} — 

| "eens 

n_aep] ~—*| real_threat|_ stance 
_ "distance 


anxiety 


EH! i 


MacCallum, R. C., Wegener, D. T., Uchino, B. N., & Fabrigar, L. R. (1993). The problem of 
equivalent models in applications of covariance structure analysis. Psychological bulletin, 114(1), 
185. 
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More equivalent models 
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ivalent models 


More equ 
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More equivalent models 
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More equivalent models 
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More equivalent models 
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What does pcalg come up with? 


Data <- read.csv(text = ' 

In-group Identification; identification; 4.56; 0.85 
Individual Deprivation;in_dep; 2.39; 0.81 

Collective Deprivation;col_dep; 3.31; 0.92 

Intergroup Anxiety;anxiety; -0.20; 0.17 

Symbolic Threat;sym_threat; 3.46; 0.76 

Realistic Threat;real_threat; 3.10; 0.88 

Personal Emotional Uncertainty;uncertainty; 2.84; 0.67 
Perceived Injustice;injustice; 2.38; 0.68 

Perceived Illegitimacy authorities;illegitemacy; 2.37; 0.02 
Perceived In-group superiority;superiority; 3.26; 0.93 
Distance to Other People;distance; 2.32; 0.66 

Societal Disconnected; disconnected; 2.79; 0.96 

Attitude towards Muslim Violence;muslim violence; 2.89; 1.06 
Own Violent Intentions;own_violent; 2.08; 0.91 

", sep = ";", header = FALSE) 

names(Data) <- c("name","var","mean", "sd") 
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library("lavaan" ) 

corMat <- getCov(' 

1 -.19 .08 -.25 .42 .07 .08 -.06 
1 .49 .36 .23 .50 .21 .50 .25 .12 
.11 .54 .62 .26 .38 .21 .31 .18 
(OT ofS Al F a2 «85 408° <~ 22-426 
64 .21 .24 .07 .39 .01 .04 .17 
27 .34 .16 .35 .19 .14 .26 .16 
.10 .08 .29 .18 .00 .30 .14 

.15 .01 .03 .23 .04 .06 

22 .17 .35 .35 .24 

.34 .08 .53 .30 

.08 .44 .39 

.24 .00 

47 


- —_ —_ —_ —_ —_ —_ —_ —_ —_— —_— —_— —_ 


, Lower=FALSE ) 
covmat <- cor2cov(corMat, Data$sd) 


-.28 .09 -.17 
wld: Bl gl 2: 209 
.09 .20 .10 
.18 .14 
-.01 


colnames(covmat) <- rownames(covmat) <- Data$var 


-.25 


-.04 


-.07 
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library ("pcalg") 

pc.fit <- pc(suffStat = list(C = covmat, n = 131), 
indepTest = gaussCItest, 
alpha=0.05, labels = as.character(Data$name) ) 
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library("qgraph" ) 
qgraph(pc.fit, vsize = 15, vsize2 = 5, shape = "ellipse", repulsion = 0.7) 


Perceived Injustice 
2 Personal Emotional Uncertainty 





Individual Deprivation 


Own Violent Intentions 





Realistic Threat) _ \ | 
KA \ 


Distance to Other People 
Perceived In-group superiority > 
Attitude towards Muslim Violence 





Symbolic Threat 





Intergroup Anxiety 


a 2 ae Societal Disconnected 
In-group Identification 
Perceived Illegitimacy authorities 
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Does it fit? 


Mod <-' 

muslim_violence ~ own_violent + distance + disconnected + superiority 
disconnected ~ identification 

superiority ~ sym_threat 

sym_threat ~ col_dep + identification 

col_dep ~ in_dep 

real_threat ~ in_dep + sym_threat 

in_dep ~~ injustice 


fit <- sem(Mod, sample.cov = covmat, sample.nobs = 131) 
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Does it fit? 


round(fitMeasures(fit)[c('chisq','df', 'pvalue','cfi','nfi','rmsea','rmsea.ci.lower','rmsea.ci.upper')], 
2) 


HH chisq df pvalue cfi nfi 
HH 80.52 39.00 0.00 0.89 0.81 
Ht rmsea rmsea.ci.lower rmsea.ci.upper 
Ht 0.09 0.06 0.12 

- Not really... 
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Directed Networks 


- Modifying a causal network until it fits can lead to over-fitting 
- Many equivalent models 


- If your goal is to test a specific hypothesis (as was the goal of 
Doosje et al.) then SEM Is good. For more exploratory research, 
different methods should be used. 


- Causal search algorithms on the other hand can lack power 
- In my experience rarely an interpretable model comes out 
- Both methods assume the true model is a DAG! 


- The arrows could be misleading, they imply a specific effect of 
doing something even though we only know what happens if we 
condition on something 


- Equivalent models are also a huge problem for network analysis 
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Background Determinants Radical Belief Violent 
Variables System Intentions 





R?=17 
Intergroup 
Anxiety 
In-group 
Identification 


Symbolic 


Individual ——— 


Deprivation i 

_R°=.18 Attitude 

Perceived towards Muslim 
In-group Violence 


Collective Realistic Superiority 
Deprivation Threat 


R°=.05 


‘ R= Distance to R2= 26 
ersona Other People j 
niki Own Violent 


Uncertainty Intentions 
R°=.05 


R2=27 Societal Dis- 
Perceived Connected 
Injustice 


Fig. 1. Final structural equation model. All paths are significant. R? = % variance explained. 


- Would, theoretically, Do( Attitude towards Muslim Violence ) really not 
influence any of the predictors? 
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Equivalent Models and Centrality 


Model A Model B Model B 
1 
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Markov Random Fields 


Markov Random Fields 


- A pairwise Markov Random Field (MRF) is an undirected network 


- Two nodes are connected if they are not independent conditional on 
all other nodes. 


More importantly, two nodes are NOT connected if they are 
independent conditioned on all nodes: 


- X; 1 X;| XO) =a) => (i,j) ¢ E 
- Anode separates two nodes if it on all paths from one node to 
another 
- No equivalent models! 
- Clear saturated model is a fully connected network 
- Naturally cyclic! 


42/110 


‘ Bseparates A and C' 
“AI C|B 
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Insomnia Concentration 


Fatigue 


Worrying and fatigue separate Insomnia and Concentration 
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Markov Random Fields in Other Sciences 


- MRF's have been used in many disciplines in Science 


Physics 

Machine Learning 
Artificial Intelligence 
Biometrics 
Economics 


Image processing 
(http://cs.stanford.edu/people/karpathy/visml/ising_example.html) 


Neural Networks (http://nxxcxx.github.io/Neural-Network/) 
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Interpreting a Markov Random Field 


The edges in a MRF can be interpreted in several ways: 


- Predictive effects 
- Arepresentation of conditional independence relationships 
- Pairwise interactions 
- Genuine symmetric relationships between nodes 
- Ising Model 
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Predictive Effects 


If this model is the generating model, does: 


' A predict B? 
- Yes! 
’ B predict A? 
- Yes! 
' A predict B just as well as B predict A? 


- Using linear or logistic regression, yes! 
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# Generate data (binary): 
A <- sample(c(0,1), 10000, replace = TRUE) 
B <- 1 * (runif(10000) < ifelse(A==1, 0.8, 0.2)) 


# Predict A from B (logistic bregression): 
AonB <- gim(A ~ B, family = "binomial") 


coef (AonB) 
## (Intercept) B 
Ht -1.443 2.873 


# Predict B from A (logistic regression): 
BonA <- glm(B ~ A, family = "binomial") 


coef (BonA) 
## (Intercept) A 
Ht -1.440 2.873 


- The logistic regression parameters are equal! 


48/110 


Predictive Effects 


‘ A predicts B and B predicts A 
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Predictive Effects 


If this model is the generating model, does: 


* A predict C’ or vise versa? 
- Yes, they are correlated 

' A predict C’ or vise versa when also taking B into account? 
- No! 


* Ina multiple (logistic) regression, C should not predict A when B is 
also taken as predictor 
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# Generate data (Gaussian): 
A <- rnorm( 10000) 

B <- A + rnorm(10000) 

C <- B + 2*rnorm(10000) 


51/110 


# Predict A from C: 
AonC <- Im(A ~ C) 


summary (AonC) 

HH 

## Call: 

#44 Im(formula = A ~ C) 

HH 

## Residuals: 

Ht Min 1Q Median 3Q Max 

## -3.521 -0.617 -0.012 0.624 3.166 

HH 

## Coefficients: 

HH Estimate Std. Error t value Pr(>|t]) 

## (Intercept) 0.01060 0.00914 1.16 Qi25 

#4 C 0.17084 0.00375 45.56 <2e-16 *** 
HH --- 

## Signif. codes: O '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
HH 


## Residual standard error: 0.914 on 9998 degrees of freedom 
## Multiple R-squared: 0.172, Adjusted R-squared: 0.172 
## F-statistic: 2.08e+03 on 1 and 9998 DF, p-value: <2e-16 
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# Predict A from B and C: 
AonBC <- Im(A ~ B + C) 


summary (AonBC) 

HH 

## Call: 

## Im(formula = A ~ B + C) 

HH 

## Residuals: 

HH Min 1Q Median 3Q Max 

## -2.9409 -0.4774 -0.0027 0.4785 2.9478 

HH 

## Coefficients: 

Ht Estimate Std. Error t value Pr(>|t|) 

## (Intercept) 0.00464 0.00711 0.65 0.51 

## B 0.49992 0.00620 80.69 <2e-16 *** 
Ht C 0.00417 0.00358 1.17 0.24 

tH s«—--- 

## Signif; codes: OQ °***" 0.001 "**" 0.01 "*" 0.05 "27 0.1°° "1 
HH 


## Residual standard error: 0.711 on 9997 degrees of freedom 
## Multiple R-squared: 0.498, Adjusted R-squared: 0.498 


## F-statistic: 4.97e+03 on 2 and 9997 DF, p-value: <2e-16 
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’ A predicts B better than C’ predicts B 
* The relationship between A and Cis mediated by B 
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Conditional Independency Relationships 


Directed Acyclic Graph 
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Directed Acyclic Graph Markov Random Field 
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Directed Acyclic Graph 
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Directed Acyclic Graph Markov Random Field 


58/110 


Directed Acyclic Graph 
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Directed Acyclic Graph Markov Random Field 
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Directed Acyclic Graph 
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Directed Acyclic Graph Markov Random Field 


- AMRF can not represent the exact implied independence relationship 
of a collider structure 


- Three edges are needed instead of two 


However, exogenous variables are commonly modeled to be 
correlated anyway 
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Directed Acyclic Graph Markov Random Field 


Easiness of Class Intelligence Easiness of Class Intelligence 
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Directed Acyclic Graph Markov Random Field 






227 


MRF's can also represent independence relationships that DAG's can 
not. 
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Data generating structure 
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Data generating structure Markov Random Field 





If we could condition on 7 
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Data generating structure Markov Random Field 





If we could not condition on 7 


Equivalent models 


- Data generated as a cluster of interacting components can fit a 
factor model perfectly! 
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The Ising Model 
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The Ising Model 


PS 
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Markov Random Fields as Generating Structure 


* In the Ising model, we could hold a magnet in some way—Do( A)— 
which can cause adjacent nodes to "flip" with the same probability if 
we conditioned on A 


- Pr(B|Do(A)) =Pr(B| A) 

- Pr(A|Do(B)) =Pr(A|B) 
- Symmetric relationship that can not be represented in a DAG 
- Real relationship that occurs in physics 
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Parameterizing Markov Random 
Fields 


The Ising Model Probability distribution 


Pr(X =a) = Sexp (s Txi + Ss" as 


<1j> 


* All X variables can typically take the values —1 and 1 


' 7; is called the threshold parameter and denotes the tendency for 
node z to be in some state 


" w4; is called the network parameter and denotes the preference for 
nodes 2 and 7 to be in the same state 
- Edge weights 


' Z is anormalizing constant (partition function) and takes the sum 
over all possible configurations of X: 


ae fe ye exp oy Te ere wij) 
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In matrix form the Ising Model probability distribution is proportional 
to: 


1 
Pr(X = x) x exp (r"a + : x' Mx) 


* g is avector of binary variables (—1 or 1) 
’ 7 is avector containing threshold parameters 


Q) is a matrix containing the network parameters and an arbitrary 
diagonal 


- Weights matrix that encodes a network 
- Wij = 0 means that there is no edge between nodes 2 and j 


- Positive weights are comparable in strength to negative 
weights 
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,T = 


—0.1 
—0.1 
—0.1 
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Pr(X = a) = Sexp (s Ti Li + ye W5jXj = 


<> 
We can compute the unnormalized probability that all nodes are 1: 
exp (—0.1 + —0.1+ —0.1+0.5+0.5) = 2.0138 


- We will call this the potential for the nodes to be in this state 


- Summing the potential of every possible state gives the normalizing 
constant Z 


- Which can then be used to compute the probabilities 
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Potential 


3.6693 
1.1052 
0.4066 
0.9048 
1.1052 
0.3329 
0.9048 
2.0138 


Z = 10.4426 


Probability 


0.3514 
0.1058 
0.0389 
0.0866 
0.1058 
0.0319 
0.0866 
0.1928 
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What if we observe Xo = 1? 
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What if we observe Xo = 1? 


%1 £2 2X3 Potential Probability 





-1 1 -l 0.4066 0.0389 
1 1 -l 0.9048 0.0866 





-1 1 1 0.9048 0.0866 


1 1 1 2.0138 0.1928 
Z = 10.4426 
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What if we observe Xo = 1? 


%1 “£2 23 Potential Probability 

=] LY «] 0.4066 0.0961 

1 1 -l 0.9048 0.2139 

-1 L 1 0.9048 0.2139 

1 1 1 2.0138 0.4761 
Z = 4.229997 


. (0.2139 + 0.4761) * (0.4761 + 0.2139) = 0.4761 
é Pr(X; =1,X3 =1| Xe =1) = Pr(X, =1| X2 =1)Pr(X3 =1| X2 = 1) 


* X, and _X3 are conditionally independent given X >! 
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Conditional Distribution for the Ising Model 


The conditional distribution for node 2 given that we observe all other 


nodes is: 
Pr (x; = 2; | x) = x () ox exp ((° = iS out “| 


j,jFt 


- This is a multiple logistic regression mode!! 


- The most common model to predict the value of one binary variable 
(dependent variable) given a set of other variables (independent 
variables) 


- The Ising model is a combination of predictive models. Edges 
represent how strongly one node predicts another 


- Apath means that the predictive strength of one node on another 
is mediated 
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Combination of Logistic models 


Os 
oS 
S 
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Ising Model as Combination of Logistic models 


(44 (4 


Pr(X, =1) «x exp(t1 + wi2%2 + wy474) 
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lsing Model as Combination of Logistic models 


(23 


Pr (Xo = 1) ox exp (TT. +wyoa1+ W323) 
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Ising Model as Combination of Logistic models 


(034 (023 


Pr (X3 — 1) x exp (73 + wo3n%9 + W344) 
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lsing Model as Combination of Logistic models 


(44 


Pr(X4=1) x exp(74 + wisr1 + w3423) 
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Continous Data 


If x is not binary but assumed Gaussian we can use a multivariate 
Gaussian distribution: 


1 1 S 
(X= 2) = exp (— 5 (ew) =e w)) 
(277) "|| 
" p iS avector that encodes the means 


* is the variance-covariance matrix 


- Now we can rearrange: 
_ f (X=) xexp(—}(@-p)'E (ep) 


_ f(X=2) «exp (- 3 (2' De —a Dt ppd tet p™S*p)) 


_ f(X=2) «exp (wie — ,2'D es) 
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Gaussian Random Field 


Reparameterizing T = To ie and Q = —S~* we obtain the following 


expression for the Multivariate Normal Distribution: 
Al a 
f(a@) xexp| 7 r+5e Oz 


Exactly the same form as the Ising Model!!! 
Except: 


- £ iS now continuous 


- The normalizing constant is different 
- The multivariate normal distribution encodes a network 


- This network is called a Gaussian Random Field (GRF), 
Concentration Graph, Gaussian Graphical Model or Partial 
Correlation Network. 
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Gaussian Random Field 


- The negative inverse covariance matrix of the multivariate normal 
distribution, also called precision matrix encodes a network 


- Through mathematical magic, standardized elements of the negative 
inverse precision matrix are equal partial correlation coefficients 
between nodes conditioned on all other nodes 


- Pig = Wij //oiitajj = Cor (Xi, X; | XO) = zi) 
- Typically these partial correlation coefficients are used in the weights 
matrix of the corresponding network structure 


- There is no edge between nodes 7 and 7 if p;; = 0 
- Which clearly corresponds to the Markov property that two nodes 


are then independent conditioned on all other nodes in the 
network 
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Conditional Distribution for the Gaussian Random 
Field 


If we observe the value of every node except node 2 we obtain the 
following (assuming all variables are standardized): 


. aia 
X;|X YON Se 25,8; 
ja i 
In which 0° is the 2th diagonal element of the precision matrix, also 
called the residual variance. 


- This is a multiple linear regression model! 


- Similar to the Ising model, the Gaussian Random Field shows how well 
nodes predict each other! 
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Combination of Linear models 


Os 
oS 
S 


91/110 


lsing Model as Combination of Linear models 


(44 (4 


—2 —2 
Xy,=71 +128, 22 +0140, 444+ €1 
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Ising Model as Combination of Linear models 


(23 


—2 —2 
X22 = T2 + W128, x2 + W2385 x3 + €2 
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lsing Model as Combination of Linear models 


(034 (023 


—2 —2 
X3 = 73 + w2383° 23 + W3483° 244+ €3 
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Ising Model as Combination of Linear models 


(44 


—2 —2 
X4 = 74 + w3404° 23 + W1404° 041 + €4 
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Markov Random Fields 


Ising Model 
- All nodes are binary 
- —lorl 


- Combination of multiple logistic regression models 
- Gaussian Random Field 
- All nodes assumed normally distributed 


- Graph structure directly related to the inverse variance- 
covariance matrix 


- Graph usually standardized to partial correlations 
- Combination of multiple linear regression models 
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Estimating Markov Random 
Fields 


Constructing Markov Random Fields 


- Because the Ising Model is a combination of multiple logistic models 
and the Gaussian Random field a combination of multiple linear 
models we can form MRF's in the following way: 


1. Predict each node from all other nodes in a multiple (logistic) 
regression 


2. This gives you two parameters for each edge in the network. For 
example: 


- The regression parameter from predicting A from B and the 
regression parameter from B to A. 


- In logistic Sea these should be about equal. In linear 
regression, they can be standardized to be about equal 


3. Take the mean of these two parameters as edge weight 


- If the regression parameter is 0 there is no edge! 
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Advanced MRF estimation 


- The methods on the previous slide work but require a lot of data! 


- To obtain a simpler interpretable and more stable network we will 
apply regularization 


- LASSO 


- Also, the assumption of multivariate normal data can be relaxed in 
MRF's 


- This will both be discussed on Thursday! 
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Examples 


Radicalization 





: In-group Identification 

: Individual Deprivation 

: Collective Deprivation 

: Intergroup Anxiety 

: Symbolic Threat 

: Realistic Threat 

: Personal Emotional Uncert 
: Perceived Injustice 

: Perceived Illegitimacy auth: 
10: Perceived In-group superi 
11: Distance to Other People 
12: Societal Disconnected 
13: Attitude towards Muslim \ 
14: Own Violent Intentions 


OONODOARWND — 
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Radicalization 


In-group Identification 
Individual Deprivation 

Collective Deprivation 
Intergroup Anxiety 

Symbolic Threat 

Realistic Threat 

Personal Emotional Uncertainty 
Perceived Injustice 

Perceived Illegitimacy authorities 
Perceived In-group superiority 
Distance to Other People 
Societal Disconnected 

Attitude towards Muslim Violence 
Own Violent Intentions 


Muslim Violence 
16.12 
20.10 
17.05 

Inf 
11.44 
14.81 
10.73 
24.83 

Inf 

3.16 
5.70 
9.15 
0.00 
ee | 


Violent Intentions 
19.89 
23.87 
20.82 

Inf 
15.21 
18.58 
14.50 
28.60 

Inf 

6.93 
8.81 
12.92 
cy | 
0.00 
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Quality of Life 


Partial Correlations (adaptive LASSO) 





@® General health 

© Physical functioning 

© Mental health 

© Role limitations physical 
© Role limitations emotional 
© Bodily pain 

© Social functioning 

© Vitality 

@® No dimension 
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Quality of Life 


Partial Correlations (adaptive LASSO) 





@® General health 

© Physical functioning 

© Mental health 

© Role limitations physical 
© Role limitations emotional 
© Bodily pain 

© Social functioning 

© Vitality 

@® No dimension 
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Quality of Life 


How true/untrue is the statement "| am as healthy as anybody | know"? 


How true/untrue is the statement "my health is excellent"? 








How true/untrue is the statement "| expect my health to get worse"? 





How true/untrue is the statement "I seem to get sick a little easier than other people"? 
In general, how would you say your health is? 





By Jolanda Kossakowski 
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Quality of Life 


Does your health now limit you in walking one block? 


Compared to one year ago, how would you rate your health in general now? 


Does your health now limit you in walking several blocks? 


Does your health now limit you in climbing one flight of stairs? 
Does your health now limit you in walking more than a mile? 


Does your health now limit you in bathing or dressing yourself? 


Does your health now limit you in climbing several flights of stairs? 


Does your health now limit you in bending, kneeling or stooping? 


Does your health now limit you in lifting or carrying groceries? 
Does your health now limit you in vigorous activities, like lifting heavy objects? 





Does your health now limit you in moderate activities, like pushing a vacuum cleaner? 





By Jolanda Kossakowski 
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Quality of Life 


Have you had difficulty performing your work or other activities? 


Were you limited in the kind of work or other activities? 


Have you accomplished less than you would like as a result of physical problems? 
Did you not do your work or other activities as carefully as usual? 


Have you cut down on the amount of time you spent on work or other activities 
as a result of physical problems? 
Have you accomplished less than you would like as a result of emotional problems? 


Have you cut down on the amount of time you spent on work or other activities 
as a result of emotional problems? 











By Jolanda Kossakowski 
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Quality of Life 


To what extent has your physical health or emotional problems 
interfered with your normal social activities? 


How much of the time has your physical health or emotional problems 
interfered with your normal social activities? 


How much did pain interfere with your normal work? 


How much bodily pain have you had during the past 4 weeks? 





By Jolanda Kossakowski 
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Quality of Life 











How much of the time during the past 4 weeks did you feel full of pep? 


How much of the time during the past 4 weeks did you have a lot of energy? 


Ly 
A 


A 
wi 
~*~ | 


Ss 


How much of the time during the past 4 weeks have you been a happy person? 


How much of the time during the past 4 weeks did you feel tired? 


How much of the time during the past 4 weeks have you been a nervous person? 
How much of the time during the past 4 weeks have you felt calm and peaceful? 


How much of the time during the past 4 weeks did you feel worn out? 


How much of the time during the past 4 weeks have you felt down in the dumps? 
How much of the time during the past 4 weeks have you felt downhearted? 


Se 
By Jolanda Kossakowski 
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Thursday: 


Estimating sparse networks 
Normality assumption 
- Another example 
- Cookbook 
- Most codes of the course! 
Be there! 
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